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Foreword 

The work described in this report is part or a continuing effort being carried 
out at UAH. The participants in this study have maintained a close liaison with 
Mr. Bob Fisher and Mr. Jerry Vanneman of NASA-MSFC, 

The purpose of this close liaison has been to discuss problems, results, and new 
techniques in order that the results of this study be as beneficial to this group 
as possible. Consequently, many informal recommendations have been made, which 
may or may not be formalized when this work is completed. 
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INTRODUCTION 


The Problem 

The problem to be studied stems from the requirement for the control system 
thrustors of a space shuttle type vehicle, necessary for exo-atmospheric flight 
control, to be buried within the heat shield of the vehicle so that interference 
heating will be reduced during flight in the atmosphere. The geometry of a typical 
installation is illustrated in Figure 1. During atmospheric flight, the non-firing 
thrustor nozzle represents a cavity in the surface of the vehicle. The interaction 
of the hypersonic external flow with the cavity produces excess local heating which 
must be ameliorated. 

Mechanical systems, such as a cover that slides over the nozzle exit, and 
active cooling systems introduce other problems besides adding undesirable weight 
and complexity to the vehicle. The simple expedient of increasing heat shield 
thickness in the region of the cavity would ultimately imply additional ablation in 
the area and probable change in nozzle size or shape from flight to flight. This may 
require additional refurbishment before each flight. 

This study, therefore, is aimed at determining the mechanisms involved in 
producing the high local heating rates in and around the thrustor nozzle in non-firing 
condition. In such a case, the nozzle is treated simply as a cavity. 

Literature Survey 

A survey of the literature available on this and related problems revealed that 
very few results applicable to this particular problem have appeared in the open litera- 
ture. Of the papers reviewed, only those referenced in this report appear to have 
direct consequence for this study. A more complete bibliography on cavity flows will 
be included in the final report for this study. 
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Experimental Program 

A test program, conducted at AEDC by personnel of General Dynamics Corp., 
has produced some results directly applicable to this study. In order to provide a 
basis for comparison or for development of empirical relations, a short study was 
made to compare the flow conditions in the wind-tunnel studies with flight condi- 
tions. Of particular concern was the question of whether the boundary layer ahead 
of the cavity was laminar or turbulent in either the wind tunnel or flight regimes. 

A map of these conditions is shown in'Figures 2 and 3. The wind tunnel conditions 

are for the von Karman Facility, Tunnel B, at AEDC^. The flight conditions are 
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fora typical reentry corridor, with no special attempt to relate the reentry trajec- 
tory to the space shuttle vehicle. 

In all of the wind tunnel tests and in the flight regime, it is expected that 
the boundary layer ahead of the nozzle cavity will be fully developed and turbulent. 
In flight, but not in the wind tunnel, the boundary layer probably should contain 
significant amounts of ablation products at the nozzle location. No attempt is made 
to account for the influence of ablation or surface heating ahead of the nozzle area, 
except to account for the total temperature profile through the boundary layer arid, 
hence, the total pressure profile. 

Since the analysis of the experimental data is being carried out at General 

Dynamics, no general discussion of the results will be made here. In so far as 

possible, comparisons of the results of this study are made with the raw data of the 
3 

experiments. These confirm the qualitative features of the flow derived from theo- 
retical analysis. 



ANALYSIS 


Fundamental Features of the Flow 

In the flight vehicle, the thrustor nozzles are to be installed in curved sur- 
faces. Presumably, the local radius of curvature of the body will be sufficiently 
large cf. the nozzle diameter so that the effects of surface curvature will be at 
most of secondary importance in this problem. The wind tunnel test program 
centered mainly on axially-symmetric nozzles set into a flat-plate model . In this 
study, only one basic configuration is considered. This is the axially-symmetric 
nozzle cavity set into a flat plate with the nozzle axis normal to the plate surface, 
as in Fig. 1 . The flow would be two-dimensional in the absence of the cavity. 

The disturbance produced by the cavity is three-dimensional, which greatly adds 
to the complexity of the problem. 

Sketches of a two-dimensional cavity flow and the central plane of the pre- 
sent problem are shown in Fig. 4. In the two-dimensional problem, the boundary 
layer separates from the upstream edge of the cavity. The dividing streamline 
leaves the separation point and intersects the downstream surface of the cavity, at 
a stagnation point. The location of the stagnation point depends upon several 
parameters such as cavity height and cavity width. But, in the steady state 2 -D 
solution, the stagnation streamline is the separation streamline. The flow within 
the cavity circulates about a vortex core, driven by the shearing stress along the 
dividing streamline. 

In the 3-D case, fluid may move in the y direction so that mass may be in- 
jected into the cavity in the central plane. In a steady state, the mass flowing 
into the cavity near the center line must flow out of the cavity at some distance 
from the central plane. The three dimensionality of the cavity clearly shows that 
stagnation should occur only in the central plane, so that the lateral flow would 
be expected to occur. The injected mass flow would increase the angular momen- 
tum of the flow in the cavity, thereby strengthening the vortex within the cavity. 
This vortex must exit the cavity with the outflow and be washed downstream. The 
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resultant vortex shape is reminiscent of the classical horseshoe vortex of a lifting 
wing of finite span. (Fig. 5) A crude model hrts been devised to indicate the in- 
fluence of this vortex on the pressure distribution in the neighborhood of the cavity. 
The model is discussed in some detail in the next section of this report including a 
comparison with experimental results. No attempt has been made to include com- 
pressibility effects, the influence of the boundary layer ahead of the cavity, or 
the interaction of the trailing vortex with the boundary layer downstream of the 
cavi ty . 

The strength of the vortex within the cavity is difficult to ascertain. Values 
may be inferred from two-dimensional results, but great care must be exercised in 
order to avoid erroneous conclusions. A discussion of two-dimensional and axially 
symmetric (annular cavity) studies is presented in a succeeding section of this re- 
port. As a result of these studies, some tentative conclusions may be drawn, par- 
ticularly in reference to future experimentation. 


Horseshoe Vortex Model 

The vortex within the cavity in the three-dimensional case must extend "to 
infinity, " close on itself, or end on a solid wall, according to a theorem due to 
Helmholtz. It has also been shown that in viscous flow, a vortex line may end on 
a solid boundary only at a stagnation point. Consequently, the vortex generated 
in the cavity in the case under consideration must have the form as described 
earlier (Fig .5). 

The horseshoe vortex must be one of the dominant features of the cavity- 
boundary layer interaction. The flow induced by the vortex would modify pressure 
distributions and velocity fields to a significant degree. Viscosity and compressi- 
bility would modify these effects, of course. 

A crude initial model was devised by considering the superposition of a uni- 
form flow on the flow due to a horseshoe vortex in steady, potential flow. Viscosity, 
compressibility and body forces are neglected. The boundary condition that there 
be no flow through the surface is imposed by the method of images. The cavity was 
not included in the geometry. 
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The velocity field induced by a horseshoe vortex lying in an xy plane at a 
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height z above the surface is given by : 
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where F is the circulation. The model is illustrated in Fig. 6. In comparison of 
this initial model with the actual flow geometry, the length of segment A-B of the 
horseshoe vortex is taken as being equal to the nozzle exit diameter. 

Since the image system utilized above guarantees that there be no flow 
through the surface z = 0, this model shows no flow into the cavity. In order to 
obtain a downwash velocity over the cavity opening, an additional line vortex 
segment is superimposed, on the flow field in the upper half plane. This new vortex 
segment is parallel to the segment AB, at a distance equal to z, below the cavity 
surface. Thus the velocity vector is obtained: 
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In this model the compressibility, the boundary conditions on the cavity walls, and 
the boundary layer upstream of the cavity have not been taken into account. The 
pressure at any point is obtained from the incompressible Bernoulli Equation. 


Vortex Strength 


Since the velocity and pressure fields in the vortex model are dependent upon 
the strength of the vortex, T , (circulation) some method for predicting the value of 
r must be specified. In the absence of any theoretical study which would specify 
F for the geometry of this problem, an appeal is made to the qualitative similarity 
between the flow in the plane of symmetry and the flow in a two-dimensional problem. 

Burggraf^ has given the following relation for the vorticity, U) , in a two- 
dimensional rectangular cavity: 


0.0615 
0 + b/a ) 


(a A ) 


2/3 


1/2 


where a is the width of the cavity, b is the depth of the cavity, L is the ’’effec- 
tive running length" to the separation point ( in our case, the distance from the 

2 

leading edge of the plate to the edge of the cavity), and F is the root-mean-square 

of the surface velocity over the perimeter of the cavity. 

2 

The values of F given by Burggraf are: 


b/a 

1/4 

1/3 

1/2 

1 

2 

3 

4 

7 

.01092 

.01732 

.03122 

.07029 

.12488 

. 15583 

. 17466 


For constant values of , L , and b, the variation of the vorticity ud with cavity 
width, a , may be obtained from 



i 
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where the subscript 1 denotes a reference condition . With a ^ taken as being 
equal to b, this variation of W with a is shown in Fig. 8. 

The heat transfer rate on the reattachmenT (downstream) wall of the cavity 
is given by Hodgson^ as: 



This variation is also shown in Fig . 8. 

Burggraf's theory is based upon the plane, inviscid, incompressible, but 
rotational, flow over a rectangular cavity. According to his theory, the heat trans- 
fer rate depends only on the cavity width-to-depth ratio. Hodgson includes com- 
pressibility . His results indicate that the heat transfer rate may vary with cavity 
width, even with b/a fixed, thus revealing a scaling effect. 

A comparison of the heat transfer rates on the reattachment wall predicted 
by these two-dimensional theories with the rates obtained in the plane of symmetry 
in the three-dimensional tests in a wind tunnel are shown in Fig. 9. Two items are 
immediately indicated by this comparison: 1) the test data is much steeper in the 

reattachment zone than that predicted by either Burggraf's or Hodgson's theory; and 
2) the tests show a decrease in reattachment heating with decreasing nozzle diameter 
as evidenced by tests with a two inch diameter nozzle, and tests with a one inch 
diameter nozzle. Hence, scaling may be a very important factor in the application 
of the wind tunnel test results to the flight vehicle. 

The vorticity predicted by Burggraf's two-dimensional theory was taken to 
predict the circulation T of the horseshoe vortex in the three-dimensional model 
discussed earlier. The pressure distribution obtained from this model are shown in 
Fig. 10. Experimental results, shown in Figs. 11 - 21, show qualitative agreement. 

Experimental Results 

Tests of the axially symmetric nozzle cavity in a flat plate at various angles 
of attack, were conducted in the AEDC-VKF Tunnel B, for various Reynolds Numbers 
and at a nominal Mach Number of 8. Wind tunnel operating conditions for the 
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various tests are listed in Table I. The flow conditions on the flat plate for the six 
cases shown in this report are listed in Table II. Isobars and isotherms in the vicinity 
of the cavity for these six cases are shown in Figs. 11 - 21 . 

The major qualitative features of the pressure field in these tests are pre - 
dicted by the horseshoe vortex model of the flow. Of particular interest in this 
respect is the high pressure zones just outside the cavity at 0 = 90°, 270 , as pre- 
dicted. The major discrepancy in trends between the experimental results and the 
horseshoe vortex model occurs just downstream of the cavity on the plane of symmetry, 
0 = 0°. Where the theory shows a decreasing pressure, the experiment shows a low 
pressure on the cavity lip and a rising pressure downstream. This discrepancy could 
be due to separation as the reattached flow tries to negotiate the sharp lip of the 
cavity exit. In any event, the agreement is remarkable in view of the crudeness of 
the theoretical model . This clearly demonstrates that the horseshoe vortex is the 
key feature of the flow field. 

One aspect of the test results is very interesting. In an effort to reduce the 
peak heating rate in the cavity, which occurs at the reattachment point, cold gas 
was bled (i .e., at low thrustor chamber pressure) through the nozzle in cases 2 and 3. 
With increasing bleed rate the stagnation point (reattachment point) moved farther 
into the cavity and the reattachment wall temperature increased. Thus the nozzle 
bleed aggravated the heating problem instead of ameliorating it. The horseshoe 
vortex model of the flow explains this phenomenon. Since for the bleed cases (as 
opposed to thrustor firing) the stagnation pressure of the bleed gas is much much less 
than the stagnation pressure of the external flow. As indicated in Fig. 22, the re- 
circulating gas is not much influenced by the bleed gas which must separate from the 
downstream (reattachment) wall just above the nozzle throat. Consequently, all of 
the bleed gas must be directed upward along the upstream wall . This injection of the 
bleed gas increases the angular momentum of the cavity flow, and therefore increases 
the vortex strength. The stronger vortex, in turn pulls in more of the external flow, 
again increasing the angular momentum in the cavity. The new stagnation streamline, 
moreover, originates higher in the upstream boundary layer so that it has a higher 
stagnation temperature than the streamlines closer to the wall, due to heat transfer 
to the plate . 



Thomke^ has studied the flow of turbulent boundary layers flowing over 
axial ly-symmetric cavities (annular grooves cut into the surface of a right-circular 
cone). The experimental tests were conducted at Mach Numbers of 2.0/ 2.5/ 3.0, 
and 4.0 and wind tunnel unit Reynold's Numbers of 0.95, 1 .02, 1 .36, 1 .43, and 
1 .30 million . A comparison of pressure distributions taken at M = 2.0, R^ ~ 0.95 
x 10^ per inch for two cavity shapes is shown in Fig. 23. 

The test results clearly show that rounding the reattachment wall produces a 
much flatter (almost constant) pressure distribution along the cavity floor. Addition- 
ally, the rounded corner results in a lower peak pressure, and the point of maximum 
pressure being moved downstream. These results would imply that rounding the 
corner would reduce heating in or near the cavity. 
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RESULTS 

At this stage of the study, conclusions may be drawn: 

1. The boundary layer-cavity interaction generates a 
horseshoe vortex reminiscent of the vortex system of 
a wing of finite span. 

2. The pressure distribution in and around the cavity is 
influenced very strongly by this horseshoe vortex. 

The influence of this vortex will extend far down- 
stream, and laterally beyond the sides of the cavity. 

3. Bleeding cold gas through the nozzle increases the 
strength of the vortex and, consequently, increases 
the peak heating rate in the nozzle. 

4. Variation of cavity geometry could tend to weaken 
the vortex, thereby relieving the heating problem. 

5. The key to passive system heat protection in the region 
of the cavity is the control of the vortex strength by 
shaping of the nozzle. 

Work is continuing to develop a means for predicting the vortex strength and thereby 
be able to calculate heating rates for the cavity flow. These means will have a 
theoretical base, but will rely heavily on experimental data. Continuation of this 
work could lead to procedures for optimizing nozzle shapes in order to find the 
optimal shape when thrustor performance and cavity heating are both considered. 
Continued support of the UAH effort is necessary in order to be able to develop 
sufficiently accurate techniques for prediction purposes. 
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Table I. Wind Tunnel Operating Conditions From 
AEDC VKF Tunnel B 


Mach 
No . 

Stagnation 

Pressure 

Stagna- 
tion 
Temp . 

Static 
Temp . 

Static 
Press . 

Freestream 
Densi ty 

Rey- 
nolds 
No . 

Free- 

stream 

Velocity 


N/m 2 

°K 

°K 

N/m 2 

kg/m 3 

10 6 /m 

m/sec 


(Ib^/in 2 ) 

ft) 

(°R) 

(lb f /in 2 ) 

(slug/ft 3 ) 

(io 6 /ft) 

(ft/sec) 

7.895 

8. 963x10'’ 
(130) 

672.2 

(1210) 

50.0 

(90) 

96.53 

(.014) 

.224xl0" 3 
(1 .399x10 j 

2.296 

(•7) 

1117 

(3665) 

7.968 

28.958x!0 5 

(420) 

716.7 

(1290) 

52.2 

(94) 

303.37 

(.044) 

.647x1 0" 3 
(4.04x10 ) 

6.562 

(2-0) 

1152 

(3780) 

8.013 

59.639xl0 5 

(865) 

, 

746.7 

(1344) 

53.9 

(97) 

599.84 

(.087) 

1.219xl0" 3 

(7.6l5xl0 _i> ) 

12.139 

(3.7) 

1180 

(3870) 
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FIG » 60 Superposition of Horseshoe Vortex and Free Stream, Using Method of 
Images to Enforce no Flow Through the Surface. 
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FIG. 7 . Horseshoe Vortex Model With Additional Vortex Segment to Induce 
Flow Through Cavity Exit Plane. 
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Fig. 10. FRIDICTED PRESSURE DISTRIBUTION BY HYDRODYNAMIC STUDY 
UNDER CASE 1 
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Fig. 18. TEMPERATURE DISTRIBUTION ( 5 R) INSIDE AND AROUND 
THE EDGE OF THE CAVITY, UNDER CASE 4 
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Fig. 20. TEMPERATURE DISTRIBUTION (°R) INSIDE AND AROUND 
THE EDGE OF THE CAVITY, UNDER CASE 5 
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Fig. 21 • PRESSURE DISTRIBUTION ( upper half, normalized 
by P-inf ) AND TEMPERATURE DISTRIBUTION '^R ( 
lower half ) UNDER CASE 6 
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